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Abstract 

We study mixed site-bond directed percolation on 2D and 3D lattices 
by using time-dependent simulations. Our results are compared with rig- 
orous bounds recently obtained by Liggett and by Katori and Tsukahara. 
The critical fractions p c site and pl ond of sites and bonds are extremely well 
approximated by a relationship reported earlier for isotropic percolation, 
(log^ ite /logp< te + \ogp c bond /\ogpl* ond = 1), where p< te and p£ ond are the 
critical fractions in pure site and bond directed percolation. 
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The mixed site-bond percolation is a natural extension of the ordinary bond 
percolation and site percolation models. The model is defined by opening sites of a 
lattice with probability p S i te and bonds with probability pbond, clusters are defined 
as combinations of open sites connected by open bonds. The system is percolating 
in a region of p s u e and pbond values in which an infinite cluster exists. 

The mixed site-bond percolation is a natural model for physical phenomena in 
randomly restricted media (accounted for by closed sites) with random iterations 
(accounted for by open and closed bonds). Introduced by Reynolds et. al. U, it 
was used in treating polymer gelation, capillary phenomena and in fracture theory 
(see and references in it) . 

Current work is concerned with the directed version of the mixed site-bond 
percolation model (Fig. 1), first introduced by Kinzel 0, |J, in which bonds can be 
open only in a certain direction, physically it can be accounted for by the existence 
of an external field, restricting interactions allowed in the system (as in the case of 
transport of particles in porous media, induced by an external electrical field), or it 
could also be applied to time-dependent phenomena, when movement in a chosen 
direction corresponds to development in time (see @, ||, 

Both pure bond and pure site directed percolation are known to belong to 
the same universality class known as DP (Directed Percolation) universality class. 
The DP universality class is known to be very robust; with a few exceptions all 
dynamical particle systems involving extinction-survival transition belong to DP 
universality class (see |7|, ||] and references in it). The mixed directed side-bond 
percolation should be expected to belong to the DP universality, although no nu- 
merical check to confirm it has ever been conducted. 

Recently Katori and Tsukahara || proposed a lower boundary for the phase 
transition line formed by p site and Pbond values, at which the percolation transition 
is taking place in the case of mixed directed percolation in 1+1 dimensions. On 
the other hand, exact upper and lower bounds for the phase transition line in 1+1 
dimensions have been recently proposed by Liggett [[10]. 

Based on a heuristic argument, Yanuka and Englman |2| proposed an inter- 
polation formula for the phase transition line for ordinary (non-directed) mixed 
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site-bond percolation in any dimension, allowing to calculate the phase diagram 
for the mixed percolation from the critical fractions p c * ite and pf ond for pure site and 
pure bond percolation, corresponding to the end-points of the diagram. Because 
of the random spread on data, the Monte Carlo results for the intermediate region 
reported by the same authors (for a number of lattices in 2D and 3D) don't allow 
to see any systematic difference from the phase transition line given by the formula, 
so that the precision of the critical line given by the formula is unclear. 

The principal aim of the current work is the evaluation of the phase transition 
line for the directed mixed site-bond percolation in 1+1 and in 2+1 dimensions 
using Monte Carlo simulations. For 1+1 dimensions we confirm the validity of 
upper and lower bounds suggested by Katori and Tsukahara M and by Liggett 



10], and, for both 1+1 and 2+1 dimensions, show that an interpolation formula, 



similar to the one suggested by Yanuka and Englman M for the ordinary mixed 
site-bond percolation, is applicable with high precision, although we are able to 
detect small systematic deviations. The time-dependent simulation approach used 
in the current work allows, simultaneously with the transition line, the evaluation 
of dynamical critical exponents for the mixed directed percolation. Our results 
imply that, not surprisingly, the mixed case belongs to the same universality class 
to which both pure bond and pure site directed percolation belong. 



To decide the critical line, we use a time-dependent simulation ||11||. Mixed 



directed percolation is regarded as a growth process, starting from a single site, the 
"origin". For 1+1 dimensions, space and time are defined as coordinates of sites 
on a square lattice perpendicular and parallel to the (1,1) diagonal, and growth 
clusters are formed by sites that can be reached from the origin moving through 
open sites and bonds in the positive time direction (see Fig. p. Similarly, for 2+1 
case we consider growth on a simple cubic lattice with time axis directed along the 
(1,1,1) diagonal. This kind of growth can be realized by updating states of sites of 
a triangular lattice (see |12[] for a discussion of the technique). We calculated the 
survival probability P(t) as the fraction of realizations containing sites connected 
to the origin at time t (realizations, surviving at time t); the average number of 
sites connected to the origin at time t, n(t), averaged over all realizations; and the 
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average square radius R 2 {t), averaged over surviving realizations. On the critical 
line, for big enough t, it is expected that these quantities are governed by power 
laws 

P(t) oc t- 5 

n(t) oc f (1) 
R 2 {t) oc t z 

The results of the time-dependent simulation are most conveniently represented 



as local slopes data (see |] | for a detailed description), for example, for the number 
of sites connected to the origin we plot T)(t) = \n[n(t)/n(t/m)]/ ln[m] against 1/t, 
taking m=7. For p site and Pbond on the critical line the true rj value is given by the 
intersection of the local slopes curve with ordinata, while in off-critical simulation 
the local slopes curves diverge for small 1/t. Survival probability P and the mean 
square radius R are treated in a similar way. Although the local slopes curves for 
R are quite insensitive to a change in p site and pbond values and thus less useful in 
determining the critical line than in the cases of n and P, we used the mean square 
radius data to check the consistency of our simulation and to estimate the z values. 

An example of local slopes curves for time dependent simulation is given in Fig. 
|2|. The data corresponds to pure site directed percolation, and allows us to support 
the critical point value provided by Onody and Neves [13[ using series technique 



as opposed to the result obtained by D.ben-Avraham et.al. [H using the transfer 
matrix technique. We only show the local slopes curves for the mean number of 
sites connected to the origin. Here, and in the rest of the current work, we take 
1000 steps and we average over 10 5 configurations. The result strongly implies 
that the the critical value obtained by Onody and Neves is correct, while the value 



provided by D.ben-Avraham et.al. |T^] seems to be off the critical point by more 
than the error margin, given in [14]. 

The overall results are given in Tab. [I] and in Tab. |2|, all Pb on d values given have 
an uncertainty of ±0.0002. The uncertainty has been determined by making sure 
that the local slopes curves for n and P, taken at pbond value, shifted from the one 
provided in Tab. [| and Tab.0 by the value of the error margin, clearly diverge as 
1/t approaches zero. 
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Repeating the argument proposed by Yanuka and Englman M in the case of 
mixed directed percolation, we have 

togPaite , l °SPbond _ , / 9 x 

where \ogpf ite and \ogpf ond are the critical fractions of bonds and sites in pure bond 
and site directed percolation. 

We compare the critical line given by (Q) with our Monte Carlo results in Fig. 
|3| and Fig. [|. The formula works surprisingly well, although if the difference 
between the interpolation formula prediction and the simulation results is plotted 
separately, as is done in Fig. [5], one can see that it can not be explained by the 
flaws of the Monte Carlo simulation. Error bars, given in Fig. |5], are calculated 
as a combination of the uncertainty of the interpolation formula prediction due 
to the uncertainty of the threshold values for pure DP cases and of the error in 
Pbond values, obtained in our simulations. The critical values for pure DP cases 
were taken as pf ite = 0.705489 ± 0.000004 [TJ and p( md = 0.644701 ± 0.000001 
in 1+1 dimensions, and for 2+1 dimensions we used the values obtained by |12j| , 
pf. te = 0.43525 ± 0.00013 and pf md = 0.38216 ± 0.00006. 

Apart from the interpolation formula, in Fig. |3] we compare the Monte Carlo 
results with exact boundaries for 1+1 mixed directed percolation, proposed by 
Katori and Tsukahara and by Liggett. Clearly, our results are in agreement with 
the bounds, presumed to be exact. 

Concerning the universality of mixed site-bond directed percolation, both in 
1+1 dimensions and in 2+1 dimensions, the dynamical exponents 5, 77 and t, ob- 
tained in our time-dependent simulations, were in agreement with the known DP 
values. For 1+1 dimensions, taking the DP exponents values given in M, the error 
margins, estimated by analyzing the spread of the local slopes data, were given 
as: 5 = 0.162 ± 0.005, 77 = 0.308 ± 0.005 and z = 1.263 ± 0.010. For 2+1 dimen- 
sions, using the DP exponents values from fli~2"| , the error margins were given as 
5 = 0.460 ± 0.020, 7] = 0.214 ± 0.025 and z = 1.134 ± 0.010. Thus, our results 
indicate that mixed directed percolation in 1+1 dimensions and in 2+1 dimensions 
belongs to the corresponding DP universality classes. 
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Concluding the results of our work, we would like to note that the interpolation 
formula suggested by Yanuka and Englman JJ] for the ordinary mixed site-bond 
percolation works surprisingly accurately in the case of the directed mixed site-bond 
percolation as well. Nonetheless, the accuracy of our Monte-Carlo results allows us 
to show a systematic difference. Our results in 1+1 dimensions are in agreement 



with exact bounds suggested by Katori and Tsukahara [[J and by Liggett [|Kj. If 
the interpolation formula could be proven to be an upper boundary for the critical 
line as our results suggest, it would be a much closer bound than the Liggett's 
bound, the best available so far. 
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Figure captions 



Mixed directed site-bond percolation in 1+1 dimensions. Percolation is al- 
lowed in directions indicated by the arrows. One can clearly see a finite 
cluster, consisting of 9 sites, connected to the origin. 

1+1 dimensions, local slopes for the mean number of sites connected to the 
origin at different values of p s u e , Pbond = 1 for all curves. Line shows the DP 
value of the exponent r], taken from [0. 

1- Psite = 0.706522 (critical value according to |T|), 2- p site = 0.7057, 3- 



Psite = 0.705489 (critical value according to [JIB]), 4- p site = 0.7053, 5- p s u e = 
0.7051. 

3. Critical line in 1+1 dimensions. Monte Carlo results are shown as points. 

1- The critical line given by Eq. |2|. 

2- The lower bound by Katori and Tsukahara ||, given as a solution of 

a 3 /3 4 - a/3 2 + 2a/3 = 1 with a = p c site and (3 = p c bond . 

3- The upper and lower bounds by Liggett [jnj, given by a = and 
a = W=fl> respectively. 

4. Critical line in 2+1 dimensions, for a simple cubic lattice. 

5. The difference between the predictions of Eq. |2| and the Monte Carlo results. 
Crosses indicate the uncertainty. 

a) 1+1 dimensions. 

b) 2+l dimensions, simple cubic lattice. 
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Table captions 



1. Critical line in 1+1 dimensions, Monte Carlo results. 



*Values for ordinary directed percolation are taken from |IJ| and [IB . 

2. Critical line in 2+1 dimensions, for a simple cubic lattice. Monte Carlo 
results. 

* Values for ordinary directed percolation are taken from [T2I . 
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Table ft] 



Pbond 


Psite 


0.644701* 


1 


0.7 


0.93585 


0.75 


0.88565 


0.8 


0.84135 


0.85 


0.80190 


0.9 


0.76645 


0.95 


0.73450 


1 


0.705489* 



Table g 

Pbond Psite 

0.38216* 1 

0.45 0.86410 

0.5 0.78725 

0.6 0.67115 

0.7 0.58770 

0.8 0.52460 

0.9 0.47510 

1 0.43525* 



